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NUMERICAL-SIMULATION METHOD FOR ROTARY METAL FORMING, 
RECORDING MEDIUM AND PROGRAM 



BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates to a method for analyzing on a computer, the plastic 
deformation process of a rotary formed body in rotary metal forming such as disk rolling, 
spinning, curling and rotary forging, as in for example the metal forming of a formed and 
fabricated materials having a ring shape in ring rolling. In particular, the invention relates 
to a technique for performing numerical simulation by modeling a rotary formed body 
using a rigid plastic finite element analysis method. 

Description of the Related Art 

Heretofore, with a background of increased speed in computer processing and 
advancement in finite element analysis methods, improvements in efficiency of process 
design is advancing due to performing virtual numerical simulation by computer before 
performing actual trials of metal working processes. 

However, with respect to the processing methods for plastic deformation by 
repeating localized contact with a rotary formed body, such as for example in the ring 
rolling process for performing metal forming of a formed and fabricated materials having a 
ring shape, multiple contact and withdrawal of the forming roll and the rotary formed body 
is repeated from the start of forming to completion of forming. Therefore there are cases 
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where in numerical simulation using the finite element analysis method, considerable 
computation time is necessary. 

For example, in the rigid plastic finite element analysis method, in a method for 
modeling the entire rotary formed body of a ring shape or the like, and moreover in a 
method for providing relatively fine finite elements in a predetermined zone inside the roll 
gap, as well as modeling the entire body (for example, N.S. KIM, S. MACHIDA and S. 
KOBAYASHI, Ring Rolling Process Simulation by the Three Dimensional Finite Element 
Method. Int. J. Mach. Tools Manufact. 30, 569 (1990) ), it is possible to reproduce highly 
accurately the complex plastic deformation in three dimensions of the rotary formed body. 
However, since the numerical simulation is also executed for the zone outside of the roll 
gap interior which has only been slightly deformed, there is the problem that considerable 
computation time is necessary. 

On the other hand, for example as with a two-dimensional cross-section 
approximation, essentially it is possible to considerably shorten the computation time with 
a method which handles the plastic deformation process in three dimensions with a two 
dimensional and an expanded two dimensional model. However, because there are many 
assumptions related to transformations for making this applicable when two 
dimensionalized, the accuracy with respect to the analysis of the complex shape is low, 
and the influence of for example the roll diameter of the forming roll and the ring radius of 
the rotary formed body can not be taken into consideration. Hence a problem arises in that 
the accuracy of the plastic strain obtained as the result of the analysis is low. 

Therefore, there is known a numerical-simulation method using a partial model for 
modeling by finite elements which subdivide only a predetermined zone for producing the 
main plastic deformation such as inside of the roll gap in the ring rolling (for example, D. 
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Y. YANG, and K. H. KIM, Rigid-Plastic Finite Element Analysis of Plane Strain Ring 
Rolling. Int. J. Mech. Sci. 30, 571 (1988) and D. Y. YANG, K. H. KIM and J. B. 
HAWKYARD, Simulation of T-Section Profile Ring Rolling by the 3-D Rigid-Plastic 
Finite Element Method. Int. J. Mech, Sci. 33, 541 (1991)). 

Incidentally, in the numerical-simulation method using the above-mentioned partial 
model, it is possible to shorten the computation time, and to reproduce highly accurately 
complex plastic deformation in three dimensions. However, if the velocity boundary 
conditions at the boundary (imaginary cutting plane) of the zone that has been modeled 
and the zone that has not been modeled are not appropriately set, a problem arises in that 
high accuracy analysis results can not be obtained. 

For example, in the numerical-simulation method according to the above- 
mentioned related art, the velocity boundary conditions are only set based on the 
assumption that the boundary cross-sections set in the rotary formed body have rigid body 
rotation in respective independent peripheral directions (circumferential directions) in only 
two cutting planes. Hence a problem arises in that, for the plastic deformation produced at 
the time of the forming process, it is not possible to reproduce for example the plastic 
deformation for where the rotary formed body is rolled out in the roll gap and elongates in 
the circumferential direction. 

Moreover, in the numerical-simulation method according to the above mentioned 
related art, since the ring rolling process being essentially a transient deformation process, 
is handled as divided steady deformation, the entire process can not be reproduced as a 
continuous phenomena. Hence a problem arises in that the accuracy of the simulation can 
not be improved. 
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DISCLOSURE OF THE INVENTION 
The present invention takes into consideration the above mentioned situation with 
the object of providing a numerical-simulation method for rotary metal forming where the 
accuracy of the numerical simulation can be improved without lengthening the 
computation time necessary, when at the time of performing numerical simulation by a 
computer for the plastic deformation of a rotary formed body at the time of a forming 
process by for example ring rolling, the boundary conditions are set based on the 
elongation deformation produced in a circumferential direction of the rotary formed body 
in boundary cross-sections set in the rotary formed body, that is, in imaginary cutting 
planes set for a partial model. Moreover an object is to provide a program for realizing on 
a computer this numerical-simulation method for rotary metal forming, and a computer 
readable recording medium recorded with this program. 

In order to achieve such an object of solving the above problems, the present 
invention: divides into a plurality of finite elements a predetermined model zone between 
two imaginary cutting planes intersecting with a circumferential direction in a rotary 
formed body rotated about an axis of rotation; expresses velocity boundary conditions of 
the imaginary cutting planes for the model zone by a function where angular velocity for 
rotation about the axis of rotation of the rotary body is made a variable; and analyses by a 
finite element method a forming process which continuously produces plastic deformation 
by localized contact with the rotary formed body. 

According to the above, with respect to transient deformation (that is, the 
deformation where the strain rate and the deformation stress change with time) where the 
diameter of the rotary formed body is increased without the axis of rotation of the rotary 
formed body being displaced, as in, for example disk rolling, cylindrical and conical 
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spinning, disk edge curling and rotary forging, at first only the partial zone producing the 
main plastic deformation is divided into finite elements. 

Then, the velocity boundary conditions are set assuming that elongation in the 
peripheral direction (for example the circumferential direction) of the rotary formed body 
is produced only inside the model zone, and the entire shape of the rotary formed body 
existing in the non-model zone which is continuous with the model zone maintains an 
approximately circular shape. Accordingly, the velocity boundary condition can be 
expressed by a function where the angular velocity for rotation about the axis of rotation of 
the rotary formed body is a variable. Hence, in comparison for example to the method for 
modeling and analyzing the entire rotary formed body, the computation time necessary can 
be considerably reduced without degrading the accuracy of the numerical simulation. 
Accordingly, it is possible to contribute to efficiency improvement of the process design of 
the metal working process, and it is possible to lower the number of trials in the actual 
metal working and to improve the stability of the quality of the product due to the metal 
working. 

In addition, in the present invention, in the case where the axis of rotation is 
displaced with addition of an increase in diameter of the rotary formed body in the forming 
process, the velocity boundary conditions are expressed by a function having three 
variables comprising; an increasing rate of the diameter of the rotary formed body, the 
angular velocity for rotation about the axis of rotation of the rotary formed body, and a 
correction value for traverse velocity of the axis of rotation. 

According to the above, in the case of a transient deformation where the diameter 
of the rotary formed body is increased and also the axis of rotation of the rotary formed 
body is displaced, as in, for example ring rolling, the velocity boundary conditions are set 
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assuming that elongation in the peripheral direction (for example the circumferential 
direction) of the rotary formed body is produced only inside the model zone, the entire 
shape of the rotary formed body existing in the non-model zone which is continuous with 
the model zone maintains an approximately circular shape, and the axis of rotation O of 
the rotary formed body is displaced for example in ring rolling, on a straight line 
connecting the axis of rotation of the main roll and the axis of rotation of the mandrel roll. 

Accordingly, the velocity boundary conditions can be expressed by a function 
having three variables comprising; the angular velocity for rotation about the axis of 
rotation of the rotary formed body, the increasing rate of the diameter of the rotary formed 
body and the correction value for traverse velocity of the axis of rotation. Hence, in 
comparison for example to the method for modeling and analyzing the entire rotary 
formed body, the computation time necessary can be considerably reduced without 
degrading the accuracy of the numerical simulation. Accordingly, it is possible to 
contribute to efficiency improvement of the process design of the metal working process, 
and it is possible to lower the number of trials in the actual metal working, and the stability 
of the quality of the product due to the metal working can be improved. 

Moreover, in the present invention, the velocity boundary conditions may be 
expressed by a function having other variables capable of conversion into the 
aforementioned variables. 

According to the above, by for example changing the coordinate system, the 
angular velocity of the rotary formed body, the correction value for the traverse velocity of 
the axis of rotation and the increasing rate of the diameter of the rotary formed body can be 
converted to other variables. 
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For example, in place of the increasing rate of the radius, the velocity boundary 
conditions can be expressed with the increasing rate of the angle with respect to the 
circumferential direction of the model zone which can be converted from the increasing 
rate of the radius, as a variable. 

In addition, for example in place of the increasing rate of the radius and the angular 
velocity of the rotary formed body, the velocity boundary condition can be expressed with 
the angular velocity of the two imaginary cutting planes being the boundary cross-sections 
of the model zone and the non-model zone which can be converted from the increasing 
rate of the radius and the angular velocity of the rotary formed body, as variables. 
Consequently, as well as being able to increase the degree of freedom of the numerical 
simulation, it is possible to improve the generality of the numerical simulation. In addition 
to this, it is possible to increase the degree of freedom of the programming, and improve 
the generality of the program. 

Moreover, in the present invention, a Euler method where the finite elements are 
spatially fixed is applied with respect to the rotation direction of the rotary formed body, 
and a Lagrange method where the finite elements are movable with the rotary formed body 
is applied with respect to a direction along the axis of rotation of the rotary formed body 
and the radial direction of the rotary formed body. 

According to the above, with respect to the rotation direction of the rotary formed 
body, by applying the Euler method which performs integration with respect to the 
variables such as the cross-section shape and the strain without moving the finite elements 
of the circumferential direction, and with respect to a direction along the axis of rotation of 
the rotary formed body and the radial direction of the rotary formed body, applying the 
Lagrange method where the finite elements are moved with deformation of the rotary 
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formed body, it is possible to follow the shape change and strain in the radial direction and 
the direction of the rotation axis with high accuracy, without altering the position of the 
nodal points dividing the finite elements of the circumferential direction. 

Furthermore, the present invention, this may involve; dividing a non-model zone 
outside of the model zone of the rotary body by elements for data storage, and of the two 
imaginary cutting planes, sequentially storing data related to the velocity boundary 
conditions output from one of the imaginary cutting planes in the elements and inputting 
from an other of the imaginary cutting planes, or computing the data for the non-model 
zone using interpolation related to the angle about the axis of rotation, based on the data at 
the one of the imaginary cutting planes and the data at the other of the imaginary cutting 
planes, and inputting from the other of the imaginary cutting planes. 

According to the above, the boundary cross-section of the model zone and the non- 
model zone in the circumferential direction, that is, information such as the shape and 
strain output from one of the imaginary cutting planes is input again to the model zone 
from the other of the imaginary cutting planes after a predetermined time delay, and by so 
doing the shape and strain of the model zone is updated. 

In this case, a Euler mesh for data retention for storing the information output from 
the model zone may be prepared in the non-model zone, and in the case where changes 
such as the deformation and strain in the non-model zone are smooth, this may be 
interpolated by a function such as for example a linear approximation. Whichever method 
is used, the accuracy when time integrating the calculated velocity field can be improved, 
and the numerical simulation can be executed with high accuracy and high speed. 



BRIEF DESCRIPTION OF THE DRAWINGS 
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FIG. 1 is a flow chart illustrating a numerical-simulation method for rotary metal 
forming according to a first embodiment of the present invention. 

FIG. 2 is a block diagram showing a structure of a simulation apparatus for 
realizing the numerical-simulation method for rotary metal forming shown in FIG. 1. 

FIG. 3 is a diagram showing a plastic deformation process of a rotary formed body 
in ring rolling. 

FIG. 4 is a diagram showing the coordinates of nodal points for a case where time 
integration in the radial direction is performed with respect to a calculated velocity field. 

FIG. 5 A is a diagram showing a stiffness matrix K obtained by introducing velocity 
boundary conditions. 

FIG. 5B is a diagram showing the stiffness matrix K obtained by expanding for 
three variables. 

FIG. 6 is a graph showing interpolation of data in a non-model zone NM. 

FIG. 7 A is a diagram showing calculation results of stress distribution in the radial 
direction in the case where a reduction in the radial direction of the rotary formed body is 
0%, being a diagram showing the results for an example according to a partial model of 
the rotary formed body. 

FIG. 7B is a diagram showing calculation results of stress distribution in the radial 
direction in the case where a reduction in the radial direction of the rotary formed body is 
0%, being a diagram showing the results for a comparative example where the entire 
rotary formed body is modeled and analyzed. 

FIG. 7C is a diagram showing a coordinate system employed in FIGS. 7 A and 7B. 

FIG. 8A is a diagram showing the calculation results of the stress distribution in 
the case where the reduction in the radial direction of the rotary formed body is 15%, being 
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a diagram showing the result for an example according to the partial model of the rotary 
formed body 3 1 . 

FIG. 8B is a diagram showing the calculation results of the stress distribution in the 
case where the reduction in the radial direction of the rotary formed body is 15%, being a 
diagram showing the results or a comparative example where the entire rotary formed 
body 3 1 is modeled and analyzed. 

FIG. 8C is a diagram showing a coordinate system employed in FIGS. 8 A and 8B. 

FIG. 9 A is a diagram showing the calculation results of the stress distribution in 
the case where the reduction in the radial direction of the rotary formed body 31 is 40%, 
being a diagram showing the result for an example according to the partial model of the 
rotary formed body 3 1 . 

FIG. 9B is a diagram showing the calculation results of the stress distribution in the 
case where the reduction in the radial direction of the rotary formed body 31 is 40%, being 
a diagram showing the results for a comparative example where the entire rotary formed 
body 3 1 is modeled and analyzed. 

FIG. 9C is a diagram showing a coordinate system employed in FIGS. 9A and 9B. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Hereunder is a description of a numerical-simulation method for rotary metal 

forming according to an embodiment of the present invention, with reference to the 

accompanying drawings. 

FIG. 1 is a flow chart illustrating a numerical-simulation method for rotary metal 

forming according to an embodiment of the present invention. FIG. 2 is a block diagram 

showing the structure of a simulation apparatus 10 for realizing the numerical-simulation 
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method for rotary metal forming shown in FIG. 1 . FIG. 3 is a diagram showing a plastic 
deformation process of a rotary formed body 3 1 in ring rolling. FIG. 4 is a diagram 
showing the coordinates of nodal points for a case where time integration in the radial 
direction is performed with respect to a calculated velocity field. FIG. 5A is a diagram 
showing a stiffness matrix K obtained by introducing velocity boundary conditions. FIG. 
5B is a diagram showing the stiffness matrix K obtained by expanding for three variables. 
FIG. 6 is a graph showing the interpolation of data in a non-model zone NM. 

The simulation apparatus 10 for realizing the numerical-simulation method for 
rotary metal forming according to the embodiment of the present invention, as shown for 
example in FIG. 2, comprises; a data input section 1 1, a velocity field calculation section 
12, a time integrating section 13, and a boundary condition modifying section 14. 
Moreover, the velocity field calculation section 12 comprises; a stiffness matrix 
constructing section 21, a boundary condition introducing section 22, a solving section 23, 
an unbalanced force calculating section 24 and a convergence judging section 25. 

The data input circuit 1 1, as mentioned below, inputs information necessary for the 
calculation of the velocity field in the model zone of the rotary formed body, such as 
information on reduction, deformation resistance and so forth. The velocity field 
calculation section 12 calculates the velocity field with respect to for example the 
deformation velocity and strain rate. The time integrating section 13 calculates the 
deformation and strain in the model zone of the rotary formed body according to the 
calculated velocity field. 

Here, in the velocity field calculation section 12, the stiffness matrix constructing 
section 21 constructs a stiffness matrix based on the information on the deformation 
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resistance, with respect to each nodal point set by the partial model of the rotary formed 
body. 

The boundary condition introducing section 22 sets the velocity boundary 
conditions (for example the displacement fixing) for the boundary cross-section of the 
modeled zone and the non-modeled zone of the rotary formed body, that is, the imaginary 
cutting plane. For example, in the case of deformation analysis, it sets a zero in the 
information for the degree of freedom of the nodal points that are not displaced in the 
stiffness matrix. 

Moreover, the boundary condition introducing section 22 is represented by three 
variables, namely: an increasing rate for the radius of the rotary formed body, a correction 
value corresponding to the traverse velocity of the axis of rotation of the rotary formed 
body, and an angular velocity corresponding to rotation about the axis of rotation of the 
rotary formed body, or a linear combination of other variables which can be converted to 
these variables. 

The solving section 23 solves linear simultaneous equations for stiffness matrices 
with intrinsic velocity boundary conditions, by for example the Cholesky's factorization 
method. 

The unbalanced force calculating section 24, as mentioned below, calculates the 
unbalanced force based on the equation for equilibrium of the force of the nodal points 
derived for example from the principle of virtual work. 

The convergence judging section 25 performs convergence judging of the 
unbalanced force, and outputs the judgment result to the time integrating section 13. Then, 
when the unbalanced force is judged to be converged, the deformation velocity and the 
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strain rate for example are calculated by time integration from the velocity vector, in the 
time integrating section 13. 

Then, the boundary condition modifying section 14 alters the velocity boundary 
conditions according to the calculated velocity vector. For example, in the deformation 
analysis, the shape of the modeled zone of the rotary formed body is changed, and for 
example the forming roll and the rotary formed body are separated at a contacting position, 
and the coefficients for constructing the stiffness matrix are changed by newly contacting 
the forming roll and the rotary formed body. 

Hereunder, is a description with reference to the figures, of the numerical- 
simulation method for rotary metal forming according to the present embodiment, for 
example the case where the ring rolling process is numerically simulated. 

Here, the ring rolling process is a process for producing a toric rotary formed body 
3 1 (refer FIG. 3) by metal forming. For example, as illustrated in FIG. 3 and FIG. 4, a 
ring-rolling mill 30 according to the present embodiment comprises; a main roll 32 
capable of rotation drive abutted against an outer peripheral surface of the approximate 
toric plate rotary formed body 31 and a rotatable mandrel roll 33 having a diameter smaller 
than for example the main roll 32, abutted against an inner peripheral surface of the rotary 
formed body 32. 

Here, the rotary formed body 3 1 is arranged so that the thick section in the radial 
direction has on either side, the main roll 32 on the outer peripheral side and the mandrel 
roll 33 on the inner peripheral side. Then, by rotationally driving the main roll 32, the 
rotary formed body 3 1 is caused to rotate, and by adjacently moving the mandrel roll 32 
towards the main roll 32 so that the mandrel roll 33 approaches the main roll 32, the thick 
section in the radial direction of the rotary formed body 31 is continuously rolled. 
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Accordingly, the rotary formed body 3 1 is such that the thickness in the radial 
direction is made thinner while maintaining the approximate toric plate shape, and the 
diameter is successively expanded. 

Next is a description with reference to the accompanying drawings, of a numerical- 
simulation method using the finite element method for the above mentioned ring rolling 
process. 

First, in step SOI shown in FIG. 1, a stiffness matrix K is generated, based on input 
data of for example, the reduction and rotation velocity of the main roll 32, the traverse 
velocity of the mandrel roll 33 and the deformation resistance of the rotary formed body 31. 

Here, in the rotary formed body 31, because a relatively large plastic deformation is 
produced in the vicinity of the roll gap 3 la between both rolls 32 and 33, and a relatively 
small plastic deformation is produced in the zone outside of the vicinity of this roll gap 
31a, partial modeling is performed by limiting the section for division into a plurality of 
finite elements for analysis to a predetermined model zone M in the vicinity of the roll gap 
31a. 

Then, the stiffness matrix K related to information such as the deformation 
resistance for example, is generated for the nodal points of each finite element. That is, 
the stiffness matrix K is used for calculating the velocity vector corresponding to the 
deformation velocity and strain rate of each nodal point. For example, the strain rate 
indicates at what degree of velocity the model zone M is distorted, and the deformation 
resistance in this case corresponds to the resistance produced when strain occurs in the 
model zone M. 

Next, in step S02, as shown in for example the following equation (1), the velocity 
boundary condition (for example the displacement fixing) is introduced to the equation of 
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the system set based on the stiffness matrix K. Moreover, in the equation (1), [v] denotes 
the velocity vector corresponding to the deformation velocity and strain rate at each nodal 
point, and [F] denotes the force vector. 

Then, as mentioned above, the stiffness matrix K does not include information 
related to for example the constraint condition of each nodal point, but is comprised of 
information related to what degree of deformation resistance each nodal point has. 

[K][v] = [F] (1) 

Here, as mentioned below, the velocity boundary conditions are set with respect to 
the predetermined model zone M in the vicinity of the roll gap 31a which is divided into 
finite elements by for example the partial model, with the two imaginary cutting planes 
3 1 A and 3 IB on either side of the model zone M as boundaries, and so as to satisfy the 
conditions where as if the overall shape of the rotary formed body 31 makes up an 
approximate toric plate shape for the reason that a continuous portion exists. 

Moreover, in FIG. 3 and FIG. 4, the rotation direction of the rotary formed body 3 1 
is a clockwise direction, and the imaginary cutting plane 3 1 A is provided on the rotation 
direction side of the two sides of the model zone M, and the imaginary cutting plane 31B 
is provided on the opposite direction side to the rotation direction. 

Then, for example in the case of deformation analysis, in the stiffness matrix K, 
zeros are set for the information of the zones that are not displaced. Accordingly, as in the 
matrix K a shown for example in FIG. 5 A, the stiffness matrix K forms a distribution A 
of a diagonal zone having a predetermined width and non-zero elements concentrated near 
the main diagonal line, thus constructing a symmetrical band matrix. 

Moreover, as mentioned below, in the stiffness matrix K to which the velocity 
boundary condition has been introduced, the velocity boundary condition is expressed by a 
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linear combination of three variables, namely; an increasing rate r v of the radius of the 
rotary formed body 31, a correction value c v for the traverse velocity o v of the axis of 
rotation O of the rotary formed body 3 1, and an angular velocity co for the rotation about 
the axis of rotation O of the rotary formed body 3 1 . 

That is, in the stiffness matrix K forming a coefficient matrix set with respect to 
each velocity component vi, ... , v n constituting a velocity vector [v] = vi, ... , v„, a zero is 
set for the coefficient corresponding to each velocity component v h ... , v n , and the three 
variables of the increasing rate r v of the radius, the correction value c v for the traverse 
velocity o v , and the angular velocity w are added to the velocity vector [v], and the 
predetermined coefficient corresponding to these three variables is added to the stiffness 
matrix K. 

Accordingly, in the matrix K /3 shown in FIG. 5B for example, the stiffness matrix 
K is constructed with a zone of three columns corresponding to the three expanded 
variables, newly reserved in the last line of the matrix K a shown in FIG. 5A. 

In this case, by utilizing for example the coefficient storage method of the skyline 
method, it is not necessary for example to expand the memory zone from the a shown in 
FIG. 5A to the /3 shown in FIG. 5B, it being acceptable to merely expand the memory 
zone of the expanded three column part. 

Then, in step S03, the solution of the non-linear simultaneous equations is obtained 
using equation (1), by for example a direct iteration method. In this case, in the situation 
where for example degrees of freedom in the three directions are provided with respect to 
each nodal point of the total number of nodal points N, this becomes a non-linear 
simultaneous equation in 3N dimensions. 
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Next, in step S04, the unbalanced force is calculated. That is, in the linear equation 
for obtaining the velocity vector [v] shown in equation (1) above, since the stiffness matrix 
[K] and the force vector [F] are expressed by functions of the velocity vector [v], firstly 
the current value [v m ] of the velocity vector [v] is calculated using the stiffness matrix [K 
(v m .i)] and the force vector [F (v m _i)] calibrated for the previous value [v m _i] of the 
velocity vector [v] in a series of processes. 

Then, based on the current value [v m ] of the calculated velocity vector [v], the 
current stiffness matrix [K (v m )] and force vector [F (v m )] are calculated, and the 
difference between [K (v m )] x [v m ] and [F (v m )] is calculated as the unbalanced force. 

Moreover, the method of solving the simultaneous equations is not limited to the 
direct iteration method above, and other methods such as for example the Newton- 
Raphson method may be used. 

Then, in step S05 it is judged if the difference between [K (v m )] x [v m ] and [F (v m )] 
has converged to be smaller than for example a predetermined threshold value, that is, it is 
judged if [K (v m )] x [v m ] and [F (v m )] retain a balance related to the force. 

In the case where the judgment result is 'NO', control proceeds to step SOI . On 
the other hand, in the case where the judgment result is 'YES', control proceeds to step 
S07. 

Then, in step S07, a time integration process is performed with respect to the 
calculated velocity field, the deformation and strain in the model zone M is calculated, and 
the velocity boundary condition is changed. 

Here, for example when the deformed shape is calculated by the time integration 
process, then as mentioned below, this integration is performed using a cylindrical 
coordinate system centered on the axis of rotation O of the rotary formed body 31 at this 
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time. Consequently, it is possible to avoid the accumulation of errors due to the 
integration calculation. 

Moreover, in the numerical integration, considering the velocity component 
corresponding to rotation in the circumferential direction of the rotary formed body 31 
becoming large, it is considered necessary to avoid that the partially modeled model zone 
M would be flowing out from within the roll gap 3 la in only the peripheral direction 
(circumferential direction). Hence the integration is performed in the circumferential 
direction using the Euler method. 

On the other hand, with respect to the cross-section shape of the rotary formed 
body 3 1 prescribed by the diametric direction (the radial direction) and the direction along 
the axis of rotation O, the integration is performed using Lagrange's method in order to 
accurately reproduce the change in shape of the cross-section form. 

Then, in step S08, it is judged if a predetermined time in the series of processes has 
lapsed. In the case where the judgment result is 'NO', the processes subsequent to step 
SOI are repeated. On the other hand, in the case where the judgment result is 'YES', the 
series of processes is terminated. 

Consequently, firstly, based on the velocity field calculated by the initial process, 
each finite element is virtually displaced towards the exit side, that is, the imaginary 
cutting plane 3 1 A. By so doing, the data (for example, the information for the deformed 
shape and strain) that is set for the imaginary cutting plane 3 1 A on the exit side before 
displacement, is output from the model zone M. On the other hand, in the imaginary 
cutting plane 3 IB on the entry side, new information that has not been set is entered into 
the model zone M before displacement. 
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That is to say, because the approximate toric plate rotary formed body 3 1 is formed 
while being rotated, it is necessary to re-enter data, such as the data that is output from the 
imaginary cutting plane 3 1 A on the exit side in the circumferential direction, into the 
model zone M from the imaginary cutting plane 3 IB on the entry side. 

Consequently, the non-model zone NM is divided by for example, the plurality of 
elements for data storage (Mesh), and the data output from the imaginary cutting plane 
31A of the model zone M is stored sequentially in each element. Accordingly, the data 
output from the imaginary cutting plane 3 1 A sequentially moves each element inside the 
non-model zone NM, and after moving all elements provided in a circumferential direction, 
is entered for a second time into the model zone M from the imaginary cutting plane 3 IB 
on the entry side. 

Moreover, in the case where the data that is output from the model zone M is 
entered a second time into the model zone M, the method is not limited to that as 
mentioned above of providing the mesh for data storage in the non-model zone NM. For 
example in the case where there is a smooth change of deformed shape and strain in the 
non-model zone NM, then as shown for example in FIG. 6, the data entered into the model 
zone M may also be calculated using the predetermined function related to the angle 6 
about the axis of rotation O, for example interpolation by linear approximation, based on 
the data for the imaginary cutting plane 3 1 A on the exit side and the imaginary cutting 
plane 3 IB on the entry side (for example, strains £ a , a b shown in FIG. 6). 

Hereunder is an explanation of the process of introducing of boundary conditions 
in the above-mentioned step S02. 

As shown in FIG. 3, since plastic deformation of the rotary formed body 3 1 is 
rendered in the roll gap 3 1 a between the main roll 32 and the mandrel roll 33, then in the 
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case where the deformation analysis is performed by numerical simulation, the partial 
model for analysis is performed by dividing only the vicinity of the roll gap 3 la into finite 
elements. 

Here, when time has lapsed from time t = to to time t' = t 0 + A t, plastic deformation 
in the rotary formed body 3 1 is produced under the following three assumptions. 

Firstly, for the first assumption, the rolling (elongation) corresponding to the 
circumferential direction of the rotary formed body 3 1 will be produced only in the inside 
of the model zone M in the vicinity of the roll gap 3 la that is partially modeled. 

In addition, for the second assumption, the non-model zone NM that has not been 
partially modeled, that is, the zone excluding the vicinity of the roll gap 31a, is 
continuously connected to the model zone M, and for the whole shape of the rotary formed 
body 31 the shape of an approximate toric plate is maintained. 

Moreover, for the third assumption, the axis of rotation O of the rotary formed 
body 31 is positioned on a straight line connecting the axis of rotation PI of the main roll 
32 and the axis of rotation P2 of the mandrel roll 33. 

Here, for example, in the ring rolling mill 30 according to the present embodiment, 
these three assumptions are reasonable assumptions because, with regards to the first 
assumption, for example the axial roll (not shown in the figure) between which the rotary 
formed body 31 is inserted in the direction along the axis of rotation O, does not perform 
significant shaping. Moreover, with regards to the second assumption, the approximate 
toric plate rotary formed body 3 1 is such that the form of the toric shape is substantially 
retained even when the diameter is successively enlarged. Furthermore, with regards to 
the third assumption, a pair of support rolls (not shown in the figure) abutted against the 
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outer peripheral surface of the rotary formed body 3 1 for restricting the displacement of 
the rotary formed body 3 1 in a predetermined direction is provided. 

Here, with respect to the radius r 0 of the circumference passing the center position 
of the thick section in the radial direction of the rotary formed body 3 1 in time t, the radius 
r 0 \ after the deformation in the time t' = t 0 + A t after the predetermined time lapse A t, is 
expressed in the following equation (2), by the increasing rate r v of the radius. 

ro'= ro + rvAt (2) 

Here, the elongation deformation in the circumferential direction of the rotary 
formed body 3 1 according to the first assumption is produced inside the model zone M. 
Moreover, the rotary formed body 3 1 according to the second assumption is such since the 
shape of an approximate toric plate is maintained, then even if the angle 0 z made by the 
model zone M about the axis of rotation O, that is the angle 0 z made by the imaginary 
cutting plane 3 1 A on the exit side and the imaginary cutting plane 3 IB on the entry side 
with respect to the axis of rotation O, is changed to angle 6 z ' in time t' after the 
deformation, as shown in the following equation (3), the length of the circumference in the 
non-model zone NM of the rotary formed body 3 1 is not changed. 

(27i-Gz)ro = (27i-e z ')ro' (3) 

Accordingly, the incremental angle A 6 z of the angle 6 z made by the model zone 
M about the axis of rotation O, is expressed in the following equation (4). 

A8„fe'-9, = r ' At(2 *- 9z > (4) 

ro + rvAt 

In addition, if the elongation deformation in circumferential direction produced in 
the rotary formed body 3 1 is divided equally with respect to the imaginary cutting plane 
3 IB on the entry side and the imaginary cutting plane 3 1 A on the exit side, and the 
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imaginary cutting plane 3 IB on the entry side and the imaginary cutting plane 3 1 A on the 
exit side are rotated at the angular velocity co about the axis of rotation O with respect to 
the circumferential direction of the rotary formed body 3 1 , then the following equation (5) 
holds between the angles with respect to the predetermined reference plane Z (shown in 
FIG. 3 and FIG. 4) that includes the axis of rotation O of the rotary formed body 31, 
namely the angle 6 a of the imaginary cutting plane 3 1 A on the exit side and the angle 9 b 
of the imaginary cutting plane 3 IB on the entry side at time t, and the angle 6 a ' and the 
angle 6 b' at the time t' after the deformation. 

e a ' = e a - — + 0>At 
2 

A0z 

8b'- 0b + — - + coAt 
2 

Here, according to the third assumption, the axis of rotation O of the rotary formed 
body 31 is displaced on the straight line connecting each axis of rotation PI and P2 of the 
main roll 32 and the mandrel roll 33. As a factor exerting influence on the displacement of 
this axis of rotation O, for example the following two factors are considered. 

The first factor is a factor for displacing the axis of rotation O accompanying the 
increase of the radius of the rotary formed body 31, in a direction away from the main roll 
32. The second factor is a factor for adjoining to the main roll 32 the circumference that 
passes through the center position corresponding to the thick section in the radial direction 
of the rotary formed body 3 1 , by adjacently moving the mandrel roll 33 towards the main 
roll 32 for which the axis of rotation PI is fixed. 

Consequently, the traverse velocity o v of the axis of rotation O is expressed in the 
following equation (6) by the radial increasing rate r v corresponding to the first factor, and 
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the correction value c v for the traverse velocity o v in a direction for where the axis of 
rotation O corresponding to the second factor approaches the main roll 32. 

Ov = Cv - rv (6) 

Moreover, as shown in FIG. 3, the traverse velocity v p i of the main roll 32 and the 
traverse velocity v p2 of the mandrel roll 33 viewed from the circumference passing through 
the center position corresponding to the thick section in the radial direction of the rotary 
formed body 31, are expressed in the following equation (7) by means of the velocity for 
moving the mandrel roll 33 towards the main roll 32, that is, the reduction velocity v p . 

Vpl = Cv ) 

(7) 

Vp2 = Vp - CvJ 

Here, with respect to the radius r n corresponding to the axis of rotation O of the 
appropriate nodal points P „ on the imaginary cutting planes 31 A and 3 IB in time t, the 
radius r n ' after the deformation in the time t' = t + At after the predetermined time lapse 
A t, is expressed in the following equation (8) by the increasing rate r v of the radius. 

rn'=r n + r v At (8) 

Furthermore, in the case where the arbitrary points (x a , y a ) of the imaginary cutting 
plane 3 1 A on the exit side and the arbitrary point (x b , y b ) of the imaginary cutting plane 
3 IB on the entry side, in the radial direction of the rotary formed body 3 1 at time t are 
moved to each point (x a \ y a ') and point (x b ', y b ') at time t' after the deformation, each 
point displacement ( A x a , A y a ) and ( A x b , A y b ) is expressed by the following equation (9). 

AXa = Xa'-Xa = rn'cOSGa'-rnCOSGa 

Ay a = ya'-ya = rn'sin Ga'-r n sin Ga + At(cv - r v ) 

■ ' \ (9) 

Axb = Xb'-Xb = rn'cosGb'-rnCosOb 

Ayb = yb'-yb = r n ' sin 0 b '-rn sin G b + At(c v - r v ) 



Each of these displacement points ( A x a , A y a ) and ( A x b , A y b ) are divided by the 
predetermined time lapse A t, and by making this time lapse A t close to zero, the velocity 
of each nodal point (u a , v a ) and (u b , v b ) is calculated, as shown in the following equation 
(10). 



v AXa ( . 2tT-9z . 
Ua = lim = COS0a + Tn Sill «a llv - Yn SU1 9a(D 

At-o At ^ 2ro 1 



sin 0 a - rn 2% 6z cos 9a - 1 |r v + r n cos GaW + c 
2ro 



Ub = lim - I cos9 

At-0 At 

, Ayb _ 



2tt -8z 

2ro 
2tt -Gz 



sin9b rv-rnSinU 



••(10) 



vb = lim— — = 1 sin9b + rn ~" "" cosGb-1 |rv + rnsin9bCD + c 
^° At V 2ro 1 



Accordingly, it is possible to express the velocity of each nodal point (u a , v a ) and 
(u b , v b ) on the imaginary cutting planes 31 A and 31 B of the model zone M, that is, the 
velocity boundary conditions, by a linear combination of coefficients comprising the 
radius r n of the imaginary cutting planes 31 A and 3 IB, and the angles 6 a and 0 b made 
by the imaginary cutting planes 31 A and 3 IB, with respect to the three variables 
comprising the increasing rate r v of the radius of the rotary formed body 31, the angular 
velocity co for the rotation about the axis of rotation O of the rotary formed body 3 1, and 
the correction value c v for the traverse velocity o v of the axis of rotation O of the rotary 
formed body 3 1 . Moreover, 8 z = 6 b - 0 a . 

Next is an explanation with reference to FIG. 4, of the process for time integration 
in step S07 mentioned above, in particular the time integration corresponding to the 
position in the radial direction of the rotary formed body 31. 

Here, in the case where the time integration of the coordinates in the radial 
direction is performed with respect to the appropriate nodal point Q 0 , the nodal point 
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adjacent to the nodal point Q G in the opposite direction of the rotating direction of the 
rotary formed body 3 1 is made nodal point Q p . 

Then, with respect to the predetermined time integration A t, the point Q e for 
reaching the position of the nodal point Q 0 after the time lapse of this time integration A t 
is positioned between the nodal point Q c and the nodal point Q p . The angle d made with 
respect to the reference plane Z about the axis of rotation O of this point Q e becomes the 
value between the angle d 0 that the nodal point Q 0 makes with respect to the reference 
plane Z and the angle 6 p that the nodal point Q p makes with respect to reference plane Z. 

Here, the angular velocity co at the point Q e is expressed in the following equation 
(1 1), by the angular velocity co 0 of the nodal point Q 0 , the angular velocity co p of the nodal 
point Q p , the velocity dr 0 /dt in the radial direction of the nodal point Q G , and the velocity 
dr p /dt in the radial of the nodal point Q p . 

»-», + ^(e-e,) (ii) 

Wo - Up 

Then, since the point Q e having this angular velocity co reaches the nodal point Q 0 
after the predetermined time integration A t, the following equation (12) holds. 

G + coAt^Uo --(12) 

Consequently, from equation (1 1) and equation (12), the following equation (13) 

holds. 



(6o-9 P ) + At(©o-G> P ) 



•(13) 
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By means of this angle 9 , the position r in the radial direction and the velocity 
dr/dt in the radial direction at the position of the point Q e is expressed by the following 
equation (14). 



0o -0 P 

£ = r = r P+ i^ ( e-e P ) 

dt Go -Gp 



••(14) 



Further, the velocity dVdt in the radial direction of the nodal point Q 0 and the 
velocity dr p /dt in the radial direction of the nodal point Q p are made the following 
equations (15). 



dr P _ 
dT~ 
dro 
dT" 



-(15) 



Accordingly, the position r 0 ' in the radial direction after the predetermined time 
integration A t with respect to the nodal point Q 0 which becomes the object of the time 
integration, is expressed by the following equation (16). 
ro'=r + rAt (16) 

As mentioned above, according to the numerical-simulation method for rotary 
metal forming of the present embodiment, the rotary formed body 31 is considered in 
order to retain the form of an approximate toric shape. By considering the behavior of the 
rotary formed body 3 1 as a whole, highly accurate numerical simulation of the ring rolling 
is possible. 

Moreover, for example compared to the method for modeling and analyzing the 
entire rotary formed body 3 1, the accuracy of the numerical simulation does not deteriorate, 
and it is possible to considerably reduce the computation time necessary. 
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Accordingly, it is possible to contribute to an improvement in efficiency of the 
process design of the ring rolling process, and to decrease the trial frequency in the actual 
ring rolling, and stability in the quality of the products made by ring rolling can be 
improved. 

Moreover, in the present embodiment as mentioned above, the velocity (u a , v a ) and 

(u b , v b ) of each nodal point on the imaginary cutting planes 3 1 A and 3 IB of the model 

zone M, that is, the velocity boundary condition, is expressed by the linear combination of 

the three variables comprising; the increasing rate r v of the radius of the rotary formed 

body 3 1 , the angular velocity co of the rotary formed body 3 1 , and the correction value c v 

for the traverse velocity o v of the axis of rotation O of the rotary formed body 31. 

However it is not limited to this, and as in a numerical-simulation method for rotary metal 

forming according to a first modified example of the present embodiment, for example, in 

place of the increasing rate r v of the radius, as shown in the following equation (17), the 

increasing rate (d Q z /dt) of the angle corresponding to the circumferential direction of the 

model zone M may be made a variable, and the increasing rate (d 9 z /dt) of this angle can 

be converted from the increasing rate r v of the radius. 

dGz ■ A0z r rv(27t-e z ) rv(27t-e z ) 

= 0z = lim = hm— i = _ ^ L (17) 

dt At-o At At-o ro + rvAt ro 

Accordingly it is possible to express the velocity boundary conditions by the 
following equation (18). 
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Vb = ro- 



Ub = ro- 



Ua = j^ro- 
Va = ro- 




sin9b - 1 rn 



2n-Q z 2 



sin0 a - 1 r n 



27T-0Z 2 



cosGb r n 



COS0a Tn 



27i-9 z 2 



2tt-Gz 2 



+ — ! 



+ 



(18) 



Moreover, as in a numerical-simulation method for rotary metal forming according 
to a second modified example of the present embodiment, for example in place of the 
increasing rate r v of the radius and the angular velocity o> in the circumferential direction of 
the rotary formed body 3 1 , as shown in the following equation (19), the angular velocity (d 
6 a /dt) and (d d D /dt) of the imaginary cutting planes 3 1 A and 3 IB may be made variables. 

Here, these angular velocities (d 6 a /dt) and (d 6 b /dt) as shown in the equation (17) 
and the following equation (20), may be converted from the increasing rate r v of the radius 
and the angular velocity co . 



Consequently it is possible to represent the velocity boundary conditions by the 
following equations (21). 





■(20) 
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Vb = -ro- 



Ub = -ro- 



Ua = ^-ro- 
( 

Va = -ro- 



sinBb - 1 ^ 

27T-0Z 



_cosOb_^ 
2tx-9z 



sinGa-1 



cosOa 
2ti-6z 



2ti-Gz 




-rnsinBb 0b 



+ r n cos0b 0b + Cv 



COS0a • 

1) 

27T-0Z 



sin0 a -K 
G 

27T-0Z 



Gb + Cv 



(21) 



As mentioned above, according to the numerical-simulation method for rotary 
metal forming according to the first modified example and the second modified example 
of the present embodiment, it is possible to express the velocity boundary conditions by 
converting the angular velocity co of the rotary formed body 3 1 , the correction value c v for 
the traverse velocity o v of the axis of rotation O, and the increasing rate r v of the radius of 
the rotary formed body 3 1 into other variables. Hence it is possible to increase the degree 
of freedom of the numerical simulation and improve the generality of the numerical 
simulation. 

However, for example as in the numerical-simulation method according to the 
partial model related to one example of the related art mentioned above (D. Y. Yang, K. H. 
KIM and J. B. Hawkyard, Simulation of T-section Profile Ring Rolling by the 3-D Rigid- 
Plastic Finite Element Method. Int. J. Mech. Sci. 33, 541 (1991)), by expressing the 
velocity boundary condition in the following equation (22), there is the case where it is not 
possible to reproduce the plastic deformation where the rotary formed body 3 1 elongates 
in the circumferential direction. On the other hand, using the proposed method, the shape 
change and strain of the rotary formed body 3 1 can be calculated with high accuracy. 
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Ua = -r n sin0a6a 

Va = ttiCOsGaGa 

. \ (22) 

Ub = -rnsinGbGb 

vb = rnCosGbGb 

Next, is an explanation with reference to FIG. 7A through FIG. 9C of the results 
for where the ring rolling process is analyzed by numerical simulation using the 
numerical-simulation method for rotary metal forming according to the above-mentioned 
present embodiment. 

FIG. 7 A to FIG. 7C are figures showing the calculation results of the stress 
distribution in the radial direction in the case where the reduction in the radial direction of 
the rotary formed body 3 1 is 0%. FIG. 7A is a figure showing the results for the 
embodiment according to the partial model of the rotary formed body 31 . FIG. 7B is a 
figure showing the result for a comparative example where the entire rotary formed body 
3 1 is modeled and analyzed. FIG. 8A to FIG. 8C are figures showing the calculation 
results of the stress distribution in the case where the reduction in the radial direction of 
the rotary formed body 3 1 is 1 5%. FIG. 8A is a figure showing the results for the 
embodiment according to the partial model of the rotary formed body 3 1 . FIG. 8B is a 
figure showing the results for a comparative example where the entire rotary formed body 
31 is modeled and analyzed. FIG. 9A to 9C are figures showing the calculation results of 
the stress distribution in the case where the reduction in the radial direction of the rotary 
formed body 3 1 is 40%. FIG. 9 A is a figure showing the results for the embodiment 
according to the partial model of the rotary formed body 3 1 . FIG. 9B is a figure showing 
the results for a comparative example where the entire rotary formed body 3 1 is modeled 
and analyzed. 
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Here, as a comparative example for the results where the entire toric plate shaped 
rotary formed body 31 is modeled and numerical simulation performed, the results for 
where numerical simulation is performed using the partial model according to the 
numerical-simulation method for rotary metal forming according to the present 
embodiment is made an example. 

Here, in the example and the comparative example, numerical simulation was 
executed with respect to the case where the toric plate shaped rotary formed body 3 1 
composed of a material with a constant deformation resistance at 100 MPa, was rolled out 
by ring rolling. 

As shown in FIG. 7 A through FIG. 9C, the stress distribution in the radial direction 
of the rotary formed body 31 is not related to the reduction, but shows a distribution 
almost the same for the example and the comparative example. Hence it can be concluded 
that even with a method using the partial model, it is possible to perform numerical 
simulation to a high accuracy commensurable with the entire model. 

However, if the computation time taken for the comparative example using the 
whole model is 1, the computation time taken for the example using the partial model is 
0.088. Hence the necessary computation time can be significantly reduced. 

Moreover, in the above mentioned present embodiment, the numerical simulation 
is performed with respect to the ring rolling process. However it is not limited to this, and 
it is also possible to apply the present invention to transient deformation where the axis of 
rotation of the rotary formed body is not displaced, as in, for example disk rolling, 
cylindrical and conical spinning, disk edge curling and rotary forging. In this case, after 
the introduction of the boundary conditions, a zero may be set for the correction value c v 
for the increasing rate r v of the radius of the rotary formed body 3 1 and the traverse 
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velocity o v of the axis of rotation O of the rotary formed body 31, and/or for the other 
variables capable of conversion into the increasing rate r v and the correction value c v . 
Consequently, as with the ring rolling, numerical simulation of rotary metal forming can 
be performed with high accuracy and high speed. 

Furthermore, in the present embodiment mentioned above, the explanation has 
been given of the process for time integration of the deformation velocity related to the 
coordinates of the nodal points. However, it is not limited to this, and it is possible to 
apply a similar integration method also with respect to obtaining other variables that are 
time integrated, such as strain by time integrating strain rate. 

In addition, it is also possible to perform temperature coupled analysis by inputting 
data for temperature dependence as a deformation resistance. 

Furthermore, the simulation apparatus 10 for realizing the numerical-simulation 
method for rotary working according to the first embodiment of the present invention may 
be realized according to specialized hardware. Moreover, this may be made constituted by 
a memory and CPU (Central Processing Unit), and a program for realizing the function of 
the simulation apparatus 10 then loaded into the memory and executed to thereby realize 
the function. 

Moreover, the numerical simulation for rotary metal forming may be performed by 
recording on a computer readable recording medium, a program for realizing the above 
mentioned numerical- simulation method for rotary metal forming according to the present 
invention, and then reading the program recorded on this recording medium into a 
computer system and executing. Further, the computer system mentioned here may 
include an OS (Operating System) and hardware such as peripheral devices. 
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In addition, a computer readable recording medium refers to a portable medium 
such as a flexible disk, a magneto-optical disk, a ROM (Read Only Memory), a CD 
(Compact Disk) - ROM, or a storage device such as a hard disk incorporated into the 
computer system. Furthermore, a computer readable storage medium also includes a 
medium which dynamically retains a program during a short time such as with a 
communication line in the case where a program is transmitted via a network such as the 
Internet or via a communication line such as a telephone line, and a medium which retains 
a program for a fixed time such as with a volatile memory inside a computer system which 
constitutes a server or a client for that case. 

Moreover, the above mentioned program may be one for realizing a part of the 
above mentioned functions. Furthermore, this may be one which can be realized in 
combination with a program where the aforementioned functions are already recorded on a 
computer system. 



